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1. Introduction 

Debye screening is a characteristic feature of any plasma. In Abelian plasmas, the electric 
field induced by a static test charge decreases exponentially with distance from the charge, 
(E(r)) oc e~"^°l'"l, with rriD denoting the Debye mass (or inverse Debye screening length). 
The Debye mass cannot be defined so easily in a non- Abelian plasma where the field strength 
is only gauge-covariant, not gauge-invariant.^ But, as discussed in Ref. [1], the Debye mass 
can be given a precise, non-perturbative definition (valid in both Abelian and non-Abelian 
theories) as the smallest inverse correlation length in symmetry channels which are odd under 
Euclidean time reflection. In the Hilbert space interpretation of a theory, Euclidean time 
reflection corresponds to the product of charge conjugation (C) and time reversal (T).^ 

In a weakly coupled plasma, the temperature T sets the scale of the energy or momentum 
of typical excitations (gluons or quarks in a QCD plasma). The Debye mass is parametrically 
smaller, mo ~ VXT, with A = g^Nc the 't Hooft coupling (defined at a scale on the order 
of T). The Debye mass mj) may be regarded as a thermally-induced effective mass for static 
fluctuations in the (Euclidean) time component of the gauge field, Aq. The equilibrium 
behavior on distance scales large compared to may be described by a three-dimensional 
effective theory which is obtained from the original four-dimensional field theory (defined 
on a thermal circle of circumference /? = l/T) by integrating out all non-static modes as 
well as the static component of Aq. What remains are just the static, spatial components 

'^For example, in the presence of a test charge (trE(r)^) is not proportional to e"^"*"'*"' at long distance. 
Rather, it falls exponentially as e^''''^^, where ^ is the longest correlation length of the system which, in weakly 
coupled plasmas, is due to static magnetic fluctuations. 

^This definition of the Debye mass is only applicable to theories which are invariant under CT, in CT 
invariant equilibrium states. Hence this definition cannot be used with non-zero chemical potentials. See 
Ref. [1] for more discussion of this point. 
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of the gauge field, whose dynamics are described by three-dimensional Yang-Mills theory^ 
with a dimensionful 't Hooft coupling A3 = AT. In non-Abelian theories, fluctuations with 
wavenumbers on this scale are intrinsically non-perturbative. Three-dimensional non-Abelian 
Yang-Mills theories are known to develop a finite correlation length equal to the inverse of 
this scale times a (non-perturbativc) pure number. Consequently, the mass gap m-gap, defined 
as the inverse of the longest correlation length, of a weakly coupled non-Abelian plasma is of 
order AT, and is parametrically smaller than the 0{\/\T) Debye mass. 

The leading weak coupling behavior of the Debye mass may be calculated from one-loop 
thermal perturbation theory. Subleading contributions to the Debye mass, and the leading 
behavior of the mass gap, depend on non-perturbative AT scale physics. In this weak coupling 
regime, inverse correlation lengths in different CT (or Euclidean time refiection) odd symmetry 
channels differ only by non-perturbative 0{\T) amounts. It is merely a convention to regard 
the Debye mass as the smallest inverse correlation length in all CT-odd channels instead of 
using the correlation length in a specific symmetry channel (such as, for example, that of the 
imaginary part of the Polyakov loop). As explained in Ref. [1], regarding the smallest CT-odd 
inverse correlation length as the Debye mass leads to a particularly simple relation between 
the next-to-leading order contribution to ttid and the expectation value of Wilson loops in 
three-dimensional Yang-Mills theory. 

The quark-gluon plasma (QGP) produced in heavy ion collision is, for much of its evo- 
lution, thought to be rather strongly coupled [2,3], and this has stimulated much theoretical 
interest in understanding the dynamics of strongly coupled plasmas. Although lattice simu- 
lations are a useful technique for extracting static equilibrium quantities, observablcs which 
are sensitive to real time evolution are not generally accessible from Euclidean lattice simu- 
lations.^ 

Gauge/string (or AdS/CFT) duality has provided new theoretical tools for studying the 
dynamics of strongly coupled gauge theories and, in particular, strongly coupled non-Abelian 
plasmas. The most accessible system is maximally supersymmetric Yang-Mills (A/" = 4 SYM) 
theory in the limit of large and large 't Hooft coupling. Although this theory is manifestly 
not QCD, at non-zero temperature AA = 4 SYM describes a non-Abelian plasma composed of 
gauge bosons and adjoint representation fermions and scalars which shares many qualitative 
similarities with hot QCD. Understanding the extent to which it is, or is not, quantitatively 
similar to a QCD plasma is of considerable current interest. 

Several quantitative comparisons of dynamic properties of A^ = 4 SYM and QCD plasmas 
in the weakly coupled regime are now available. These include the shear viscosity [11], the 
heavy quark energy loss {dE/dx) [12], and the photo-emission spectrum [13]. An interesting 
common feature has emerged from these weak-coupling comparisons: properties of AA = 4 
SYM agree surprisingly well with those of QCD provided one compares the two theories not 

^Up to higher dimensional irrelevant operators. Scalar fields, if present, generically receive thermally 
induced masses comparable to the Debye mass, and may also be integrated out. 

''However, one can attempt to constrain suitable models of spectral functions by applying maximal entropy 
fitting methods to numerical data for Euclidean correlators. See, for example, Refs. [4-10] . 
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at coinciding values of the gauge coupling, but instead compares at coinciding values of the 
Debye mass (or other closely related thermal scales) . 

One would like to understand the extent to which properties of strongly coupled A/" = 4 
SYM plasma mimic those of a QCD plasma at temperatures of perhaps 1.5-3 where the 
QCD plasma is thought to be relatively strongly coupled. In detail, this will undoubtedly 
depend on what observables (or ratios of obscrvablcs) arc considered. 

As one step in this direction, our goal in this paper is to discuss what is known about 
the Debye mass and related phenomena associated with screening in strongly coupled A/" = 4 
supersymmetric Yang-Mills theory (via AdS/CFT duality) and to compare with available 
data from lattice QCD. 

Polyakov loop (or Wilson line) correlators provide the simplest probe of screening in hot 
gauge theories. A number of previous papers, starting with Rcfs. [14,15], have used AdS/CFT 
duality to examine the behavior of Polyakov loop correlators in strongly coupled M = 4 SYM 
at non-zero temperature. Confusing physical interpretations have appeared in several papers 
(such as the suggestion that the connected part of the correlator drops abruptly to zero at 
a finite separation [16]), and none of these papers specifically address the non-pcrturbative 
definition of the Debye mass mentioned above. Therefore, we revisit the analysis of Polyakov 
loop (and related) correlators via AdS/CFT duality. In particular, wc emphasize that the 
connected parts of these correlators must be non-vanishing at all separations and do have 
exponential long distance behavior.^ 

We summarize what is known about the resulting correlation lengths in different symme- 
try channels, and compare with available data for corresponding correlation lengths obtained 
from lattice simulations in QCD at temperatures above, but near, the confinement /decon- 
finement transition temperature. An appendix discusses the transition from weak to strong 
coupling in more detail, including partial results on the first subleading 0(A~^/^) strong 
coupling corrections to the Debye mass. 

It should be noted that some authors (for example, Ref. [18, 19]) refer to "screening" 
when discussing the behavior of the static quark potential at no n- asymptotic distances, and 
interpret the phrase "screening length" to mean some (loosely defined) notion of when the 
potential deviates significantly from its zero temperature form. This is not what we mean 
by screening. We will always use "Debye screening" in the conventional sense, in which a 
screening length characterizes the asymptotic behavior of suitable correlators. 

2. Polyeikov loop correlator from AdS/CFT 

For large and large 't Hooft coupling A, the expectation values and correlation functions of 
Wilson loop operators in the AA = 4 theory can be calculated using the AdS/CFT correspon- 
dence. The path integral in the five dimensional gravity description in this limit is typically 

closely related discussion with similar conclusions for the case of topologically trivial Wilson loop 
correlators at T = can be found in Ref. [17]. 
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b) 



c) 



Figure 1: String configurations contributing to the Polyakov loop correlator. The upper lines rep- 
resent the boundary of the geometry whereas the dashed lower lines denote the location of the black 
hole horizon. The periodic temporal direction is suppressed. Configuration a) depicts a disconnected 
worldsheet, with two independent strings running straight from the boundary to the horizon. This 
gives the disc;onncc;tcd part of the Polyakov loop correlator. Configuration b) is a smooth connected 
worldsheet which generates the dominant contribution to the connected correlator for sufficiently small 
separations. Configuration c) is a connected contribution obtained by joining two otherwise discon- 
nected horizon-crossing string worldsheets by a small tube representing the graviton propagator. This 
gives the dominant contribution to the connected correlator for sufficiently large separations. 

dominated by a semi-classical open string worldsheet whose boundary coincides with the 
loop [20,21]. In this framework, the correlation function of two Polyakov loops. 



was first obtained in Refs. [14, 15]. From this, one may extract the static quark-antiquark 
potential (or more properly, free energy) at finite temperature in the usual fashion. 



In a non-confining phase (such as A/' = 4 SYM at any T > 0), the correlator (2.1) receives 
both connected and disconnected contributions. The disconnected contribution, depicted in 
Figure la, corresponds to two straight strings stretching from the boundary to the horizon. 
For sufficiently small separations |x|, the connected part of the correlator is dominated by a 
smooth worldsheet of a string connecting the two Polyakov loops. This is illustrated in Fig. lb. 
The regulated^ action of this worldsheet is finite and leads to a quark-antiquark free energy 

®To regulate the infinite worldsheet area, Refs. [14, 15] subtracted the on-shell action of two straight strings 

stretching from the horizon to the boundary. This corresponds to rcnormalizing Polyakov loops in such a 
way that their thermal expectation value is precisely A^c- (We define the Polyakov loop P(x) as a trace 
of the holonomy, not the trace divided by A'c, so its expectation value in a deconfined phase is of order 
A^c-) Strictly speaking, this prescription is not quite right, as it corresponds to a temperature dependent 
renormalization condition. Temperature independent renormalization (such as holographic renormalization 
[22]), would produce a result for the correlator (2.1) which differs by an overall multiplicative factor. Such a 
finite multiplicative renormalization is irrelevant for our purposes, but its presence would complicate portions 
of the following discussion. Therefore, we will also adopt the simple subtraction scheme used in Ref. [14, 15]. 



C(x)^(P*(x)P(0)), 



(2.1) 




(2.2) 



with attractive Coulombic short distance behavior. This connected worldsheet configuration 
ceases to exist (as a real solution) beyond a critical separation |x| > rmax = 0.277/r. Even 
before that, at a separation r* = 0.240/r, the regulated action of the smooth connected 
worldsheet (which is negative for r < r*) crosses zero and becomes positive for r* < |x| < rmax- 
As we discuss below, the connected Polyakov loop correlator, 

Cconn(x)^(P*(x)nO))-|(7')^ (2.3) 

must be a convex function which decreases monotonically with increasing separation. It can- 
not possibly drop to zero abruptly at some critical separation. Since the smooth worldsheet 
configuration of Fig. lb ceases to exist when |x| > rmax, there must be some other configura- 
tion which yields the dominant contribution to the connected correlator in this regime. 

To identify the missing contribution, it is helpful to recall the appropriate Nc counting 
(as was also explained in the closely related case of correlation functions of two topologically 
trivial Wilson loops in Ref. [17]). The disconnected part of the correlator C(x) scales as 
when ATp — > oo. Intuitively this factor of can be seen in the supergravity picture by noting 
that the endpoints of the string at the horizon have a free Chan-Paton label indicating which 
D3 branc the string ends on inside the horizon. Hence there are really configurations of 
two independent strings, as depicted in Fig. la.'^ 

Large Nc factorization implies that the disconnected part of the correlator C(x) is the 
only order contribution; the connected part is 0{Nc), or smaller by a factor of 1/iV^. Any 
connected worldsheet, such as the smooth connected worldsheet found in Refs. [14, 15], and 
sketched in Fig. lb, does indeed produce an order one contribution to the correlator. (Hence, 
it makes no sense to compare the free energies of the string configurations shown in Fig. la 
and lb, as they don't enter at the same order in A'^c-) 

What then is the correct string configuration which yields the leading connected contri- 
bution to the Polyakov loop correlator at large separation? The answer must be a connected 
worldsheet. Since, for any jxj > rmax, no such configuration exists solving the equations fol- 
lowing from the classical Nambu-Goto action, it must be a configuration with large worldsheet 
curvatures so that quantum fluctuations of the worldsheet become important and stabilize 
the configuration. When |x| ~ rmax, it is difficult to find the resulting configuration explicitly. 
But when jxj ^ rmax, the relevant worldsheet will approach a configuration of two straight 
strings hanging straight down toward the horizon, connected by a very thin tube, as is illus- 
trated in Fig Ic. This corresponds to the exchange of the lightest supergravity mode coupling 
to both strings. 

The emission and absorption vertices for supergravity modes each come with a factor of 
the square root of Newton's constant, \/Gjv ~ l/-^c, so the connected diagram with the super- 
gravity exchange does indeed contribute at the same order as does the smooth connected 

'^Equivalently, in the Euclidean description of the AdS-Schwarschild background, where there is no horizon, 
each connected string worldsheet with k holes and I handles contributes at order (A^c)^"*'"^'. So the two 
separate worldshcets, each with a single boundary hole and no handles, in the disconnected contribution yield 
an overall 0{N'^) result, while a connected worldsheet with two boundary holes is 0{N^). 
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worldsheet in Fig. lb. For large separations, it is this "connected by graviton exchange" 
diagram that gives the leading connected contribution to the Polyakov loop correlator. The 
result will scale as e~"*l''l, where m is the mass of the lightest supergravity mode that can be 
sourced by the strings.^ 

The real part of the Polyakov loop is even under CT (or Euclidean time reflection), while 
the imaginary part is C7"-odd.^ As in the field theory discussion of Ref. [1], when considering 
the correlator of just the imaginary part of Polyakov loops, 

C(x) = (ImP(x) lmV{0)) , (2.4) 

the symmetries of the operator ImV will restrict which modes can be exchanged. At large 
separations, the correlator C(x) will be dominated by exchange of the lightest mode in the 
supergravity multiplct which can be sourced by the imaginary part of the Polyakov loop. Let 
m denote the mass of this mode. Which CT-odd modes can contribute, and the resulting value 
of m will be discussed below. Suitably distorted Polyakov loops can couple to all symmetry 
channels. Therefore, what is clear on general grounds is that the Debye mass (defined as the 
smallest CT-odd inverse correlation length) is determined by the lightest CT-odd supergravity 
mode, while the ImV correlator has an exponential tail, (7(x) ~ e^™l^l, with a mass m > ttid- 
To recap, the lightest mass of all supergravity modes determines the mass gap rrigap (or 
inverse of the longest correlation length) of the M = 4 SYM plasma, while it is the mass of 
the lightest CT-odd supergravity mode which determines the Debye mass mo- We will see 
shortly that the lightest supergravity mode is CT even, and hence the Debye mass mo is 
different (and larger) than the mass gap rrigap. Both of these masses will be finite in the 
A ^ oo limit, and of order T. As discussed in Ref. [1], one can also consider a more refined 
classification which distinguishes correlation lengths in all possible symmetry channels. In 
the strongly coupled M = 4 SYM plasma, these correlation lengths will be determined by the 
lightest masses of supergravity modes of a given symmetry. This will be discussed explicitly 
below. 

The full Polyakov loop correlator (in the large A^c and large A regime) is the sum of the 
disconnected, smooth worldsheet, and high curvature horizon-crossing worldsheet contribu- 

*This is the same point which was made in Ref. [17], where it was argued that Wilson loop correlators at 
large separation are dominated by the lightest supergravity mode. There was also an early attempt [23] along 
these lines to identify the Debye mass for the J\f = 4 plasma. However, these authors considered the lightest 
mode in the spectrum which gives the mass gap instead of the Debye mass. In addition, at the time that paper 
was written the full supergravity fluctuation spectrum had not yet been worked out and it was incorrectly 
believed that the lightest mode would come from the dilaton. Only later was it realized that the lightest mode 
is actually part of the graviton. 

®We have glossed over an irrelevant subtlety. Strictly speaking, the string configurations we have been 
discussing represent contributions to the correlators of Maldacena-Polyakov loops, which differ from ordinary 
Polyakov loops by the inclusion of a linear combination of scalar fields in the exponent [20] . However, one may 
choose this linear combination of scalars to be even under both C and T, so that the imaginary part of the 
loop corresponds to a CT-odd operator. It should also be noted that we are assuming that one is considering 
the pure phase of the dcconfined SU{Nc) N = 4 SYM plasma in which the expectation value of the Polyakov 
loop is real and positive. See Ref. [1] for more discussion of the significance of this point. 
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tions depicted in Fig. 1, 



(2.5) 



With the renormalization prescription described in footnote 6, the disconnected contribution 
has a fixed value, 

Ca = . (2.6) 

For |x| < r*, the dominant connected contribution comes from the smooth worldsheet and 
yields 

C6(x) = 6-^-4(1x1) ^ (2.7) 

where ^(|x|) is the regulated area of the smooth worldsheet configuration of Refs. [14, 15], 
measured in units of the AdS curvature radius. The prefactor of VX/ (27r) is the tension of the 
string in AdS units. In this regime, the regulated area is negative, so the connected correlator 
is exponentially large compared to unity. 

For |x| > r*, the regulated area of the smooth worldsheet is positive and its contribution 
is exponentially suppressed (in y/X). The dominant contribution to the connected correlator, 
in this regime, comes from the horizon-crossing string worldsheets involving thin tubes with 
high curvature at their ends. (A transition region whose width scales as A~^/^ connects these 
two regimes. This is discussed below.) For |x| S> r*, this contribution to the correlator falls 
exponentially with increasing separation, 

Cc(x)~ce-'"«^pW, (2.8) 

where mgap is the lightest mass of all supergravity modes. The CT-odd correlator C(x) will 
have faster e~™l^l fall-off for large separations, while for small separations |x| < r* it differs 
from Cconn(x)/2 only by exponentially small (in a/A) contributions (arising from horizon- 
crossing worldsheets connected by a thin tube). 

Although it is not strictly necessary for the subsequent portions of this paper, it is 
instructive to make a plot of the Polyakov loop correlator, and the resulting quark-antiquark 
potential (2.2), for various values of A. Fig. 2 shows the connected part of the correlator for 
A = 10, 10^, 10^ and 10^, while Fig. 3 shows the resulting static quark-antiquark potential, 
given by^°'^^ 

Vqg{^) = -p-' ln[C„ + a(x) + Ce(x)] , (2.9) 

^"Note that, because the static quark-antiquark potential involves the logarithm of the full correlator and 
not just its connected part, Vqq(x) is not simply proportional to the regulated area .4(x) when r < r*. Only 
when .A(|x|) <C —4-k(IiiNc)/V^ does the static potential reduce to P~^^A{\x.\). 

These plots arc qualitatively, and even scmi-quantitatively, correct. But they arc not (and cannot be) 
exact for several reasons. First, only the asymptotic form (2.8) of Cc(x) is known. As we discuss in the next 
section, the value of the mass gap can be extracted from known results for linearized gravitational fluctuations 
in the AdS-Schwarschild background (and one finds mgap — 2.3361 ttT). But neither the 0(1) amplitude c, 
nor corrections to the asymptotic form, are known. To make a qualitatively correct figure, we have assumed 
that Cc(x) is a pure exponential at all distances, with a value of unity at rmax- Second, the regulated smooth 
worldsheet area .4(x) is ordjr defined (as a real function) for |x| < r,nax- For r, < |x| < rmax, the contribution 
Cj,(x) is a nonperturbative, exponentially small correction to the correlator. Nevertheless, when plotting the 
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Figure 2: The connected part of the Polyakov loop correlator Cconn(x) plotted as a function of the 
separation |x|, for A = 10 (blue, dotted), 10^ (green, short dash), 10'^ (red, long dash), and 10'* (black, 
solid). As A ^ oo, the correlator develops a kink at |x| = = 0.240/T. At separations larger than 
r.^,, the exponential decay of the correlator is nearly independent of A. 

for the same set of A values. As these plots clearly illustrate, as A increases the correlator, 
and the static potential, develop a kink at r = r,,,, with exponential sensitivity to A at smaller 
separations and almost no sensitivity to A at larger separations. 

If one considers the scaling with A of the static quark-antiquark potential, or the force 
Fqq(x) = —Wqq{x) between a static quark and antiquark, then (outside the narrow transition 
region) for r < r^, the force is 0(\/A) while for r > r* the force is 0{1). Focusing only on 
this comparison might lead one to think that the r > r^, contribution should be regarded as 
a subleading 1 /A-suppressed correction to the leading large A contribution from the smooth 
worldsheet. This perspective is, however, highly misleading. The smooth worldsheet (Fig. lb) 
and the high-curvature horizon-crossing worldsheet (Fig. Ic) represent different stationary 
points of the worldsheet effective action. Except within the narrow transition region, one 
configuration gives the leading large A contribution to the connected correlator Cconn(x), while 
the other gives an exponentially suppressed relative contribution. But which configuration 
is leading and which is exponentially suppressed switches at the cross-over separation r^. 

correlator for any fixed finite value of A, one must decide how to treat this contribution when |x| passes 
''max- Our lack of knowledge about the right answer (which might be computable using suitable analytic 
continuation in the string worldsheet functional integral) amounts to a non-perturbative uncertainty which, 
formally, is much less important than the (unknown) corrections to the dominant contribution Cc(x) which 
are suppressed by inverse powers of A. But if one regards C6(x) as abruptly dropping to zero when |x| exceeds 
Tmax then, for any finite A, this introduces a discontinuity in the correlator (and a delta-function spike in the 
force Fqq{-x) = —\'Vqq{x)) which is obviously unphysical. To avoid such spurious artifacts (and solely for the 
purpose of making a qualitatively correct plot), we made an ad-hoc but smooth extrapolation of ^(x) beyond 
'"max- Specifically, for \x\ > rmax we used an exponential function of the form C6(x) = a + /3e~^'^' with a, (5 
and 7 determined by matching the value and first two radial derivatives of C(,(x) at rmax- 
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Figure 3: The static quark- antiquark potential (or free-energy) Vgg(x) plotted as a function of the 
separation |x|, for iVc = 3 and A = 10 (blue, dotted), 10^ (green, short dash), 10'^ (red, long dash), 
and 10"* (black, solid). As A ^ oo, the static potential develops a kink at |x| = r» = 0.240/r. At 
separations larger than r*, the potential continues to rise in a manner nearly independent of A, and 
asymptotically approaches — TlniV^. 

Obviously, the contribution of Fig. Ic has nothing to do with 1 /A corrections which will arise 
from small fluctuations around the smooth string worldsheet of Fig. lb. 

As a final point of this section, let us clarify the nature of the cross-over near r=K . In any 
thermal field theory, consider the Euclidean two point correlator. 



of any operator W, where x is a spatial separation between the two operators. Finite tem- 
perature implies that one direction — which one normally regards as (Euclidean) time — 
is compactified. But one may equally well choose to regard the compactified direction as 
a spatial direction, and view the separation x as defining the Euclidean time direction. A 
Hilbert space interpretation based on this definition of time immediately yields the spectral 
representation. 



for the correlator. Here {en} are the excitation energies of eigenstates {\n)} of the Hamiltonian 
(for the spatially compactified system) and Cn = {n\W\0). Since every term in the sum is 
positive, this representation shows that the correlator ^(x), as a function of separation |x|, 
is positive, monotonically decreasing, and convex. Moreover, if the connected part of the 
correlator, ^(x) — ^(oo), is non-zero at some separation x, then it must be non-zero (and 
positive) at all x. In other words, a connected correlator cannot be non-zero for some range 
of separations, and then vanish identically beyond a critical separation. 



g(x) = {w\x)w{o)) 



(2.10) 





n 
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These general principles are completely consistent with the above description of the 
Polyakov loop correlator C(x). But the change in behavior near r* implies that two very 
different energy scales contribute to the spectral representation for C(x). Specifically, there 
are states with 0{T) energies starting at rrigap, and in addition states with 0{^/\T) energies 
starting at a threshold of 7?-\/A/(27r), where r/ = 20.7 T is the slope with which the regulated 
area of the smooth worldshcct crosses zero at |x| = r*, v4(|x|) ~ ?7(|x|— r*). For any finite 
value of the 't Hooft coupling A, the correlator will be a smooth convex function of separation. 
The low energy states completely dominate when |x| is big compared to r* while the high 
energy states dominate for |x| well below r*. High and low energy states make comparable 
contributions only in a narrow cross-over region around r*, whose relative width Sr/r^, scales 
as A~^/^. In the A ^ oo limit, this width shrinks to zero and the correlator develops a kink, 
but remains positive, monotonic, and convex. 

3. Correlation lengths from supergravity modes 

Fortunately the analysis of the complete TIB supergravity spectrum, in i?-symmetry singlet 
channels, on the AdS black hole background has already been performed in Ref. [24] (fol- 
lowing earlier work on portions of the spectrum in Ref. [25]). In those papers the spectrum 
was presented as glueball masses of a three-dimensional confining gauge theory (containing 
Kaluza- Klein towers whose spacing is comparable to the mass of the lightest glueballs). We 
find it much more natural to interpret the supergravity spectrum as giving us information 
about the finite temperature J\f = 4 SYM plasma. As just discussed, the supergravity spec- 
trum directly yields the spatial correlation lengths, in different symmetry channels, of the hot 
plasma. 

For the reader's convenience, we reproduce in Table 1 the results for the supergravity 
modes from Ref. [24]. The relevant supergravity fields are the metric G^jy, the dilaton-axion 
pair e~'^ -|- ia, and the NSNS and RR two forms -B^jy and C,j,y. For each one of these fields, 
the mass"*^^ of the lightest mode with a given polarization is listed, in units of vrT, together 
with its transformation properties J^^'^ under (3-|-l)-dimensional Poincarc symmetry. In 
this table (as in Ref. [24]), the longitudinal direction defined by the spatial wave vector k of 
the mode is taken to point in the direction. For our application to two-point correlators, 
this corresponds to the direction of the spatial separation x of the operators. The indices 
i and j denote transverse spatial directions (corresponding to xi and X2), and a refers to 
the directions along the internal S"^ . The table also lists the quantum numbers * for the 
0(2) X Z2 Euclidean symmetry in the x, y, t plane transverse to the longitudinal direction, 
plus charge conjugation, which is the relevant symmetry group for a transfer matrix in the 
x^ direction. Here, ^ is the two-dimensional angular momentum, Rt is the eigenvalue corre- 

^^This terminology is a bit sloppy. The mass mo is the smallest magnitude of an imaginary wavevector k for 
which the linearized supergravity equations for the various fields (in the Euclidean AdS black hole geometry) 

admit solutions whoso only spacotimc dependence is of the form exp(fi; ■ x)/(tt), where u is the radial AdS 
coordinate. This is precisely what is needed to determine the correlation lengths of Euclidean correlators. 
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Table 1: Spectrum of i?-singlet IIB supergravity modes in units of ttT, together with SYM operators 
dual to the indicated mode (see text for details). J^'^'^ are quantum numbers for (3+l)-dimensional 
Poincare symmetry, while d%^* are quantum numbers for the 0(2) x Z2 Euclidean symmetry in the 
a;yt-plane transverse to the X3 direction (plus charge conjugation). 

sponding to the reflection t — > —t, Ry is the eigenvalue for the reflection y — > —y, and C, as 
usual, is the eigenvalue for charge conjugation. (Note that only one-dimensional irreducible 
representations with 3 = need to be labeled with Ry. Two dimensional representations 
with 5 7^ contain components with both signs oi Ry.) As mentioned in the Introduction, 
Euclidean time reflection is the product of time reversal and charge conjugation, so Rt = CT. 

The field Bfj_u couples to the current Jp'^ of the fundamental string which is electrically 
charged, while the 0^,^ field couples to the magnetically charged D-string current. Hence 
these fields are both charge conjugation odd, but have opposite parity and time reversal 
assignments. 

Using the standard dictionary of the AdS/CFT correspondence, all operators dual to the 
above supergravity modes can be identified. In particular, the operator dual to the graviton 
corresponds to the energy-momentum tensor T^^, of AA = 4 super Yang-Mills theory. ""^^ The 
axion is dual to the Pontryagin density | tiF'^^F^i, = ii E ■ B, whereas the dilaton is dual to 
the J\f = 4: Lagrange density ^ = \ ^^F^'^F^^ Since the C and T quantum numbers of 

E and B are — h and , respectively, one finds H — for E ■ B, and -|— |- for £., in agreement 

with the listed assignments for the axion and dilaton modes. The operator dual to the lowest 
mode of involves the trace of four powers of the field strength; a more explicit form may 

^^Thc analysis of Rcf. [24] used a gauge ehoice in whieh the independent modes of the metric and two-form 
fluctuations correspond to polarizations without any indices in the longitudinal {xs) or AdS radial directions. 
Hence the absence of these field components in the table. Because T^i, is conserved, the corresponding "missing" 
components of the energy-momentum tensor (Tos, Tis and T33) are not independent from the components shown 
in the table. 
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be found in Ref. [26]. The two form fields B^i, and C^i, couple to the SYM antisymmetric 
tensor operator O^i^ [27], which is the J\f = A completion of the QCD operator 



O^i, = tr 



F^^aF^f^F.p + I F^^F^^^F^f^ . (3.1 



Five dimensional antisymmetric tensor fields have in general six real independent on-shell 
degrees of freedom. But in type IIB supergravity on AdSs x S^, there is an extra constraint 

linearly relating Bf^^ and Cfj^u [24] . Hence the total on-shell degrees of freedom in these two 
fields arc reduced to only six. Among these, Bij is dual to Oij and Cij to e.jj03 Co3- One can 
also find that BiQ is dual to Oio while Cjo is dual to ejosj O^j. Since B^y + iC^p or e~'^ + ia are 
multiplets of the 5.^(2, R) symmetry of the tree-level IIB supergravity, their spectrum should 
be degenerate. Indeed this can be checked from Table 1. Note that the different correlation 
lengths for 0\2 and O^j reflect our choice of X3 as the longitudinal direction. 

Examining Table 1, one sees that the true mass gap of the A/" = 4 SYM plasma (in the 
A ^ 00 and Ac ^ 00 limits) arises from the time-time component of the graviton. This 
belongs to the CT even sector, as asserted earlier. The SYM mass gap is thus 

mgap = 2.3361 (ttT) , (3.2) 

and this characterizes the longest correlations in this non-Abelian plasma. On the other hand, 
the lowest mass in a CT-odd channel comes from the response of the axion field a. Therefore 
the Debye mass for the A/" = 4 SYM plasma in the strong coupling limit is given by 

mo = 3.4041 (ttT) . (3.3) 

The imaginary part of the Polyakov loop has J^^^ = 0^ or 0^^* = 0^^ and, in particular, 
is C-odd. Hence exchange of the C-even (but T-odd) axion field cannot contribute to the ImV 
correlator. Moreover, as Tabic 1 shows, there is no supergravity field with 0"' ^ quantum 
numbers. However, the three-form field strength ((iC)i23 does have precisely this symmetry, 
and its spectrum is the same as that of C12 for our choice of longitudinal direction X3. 
Consequently, exchange of the two- form C can contribute to the ImT^ correlator, implying 
that^^ 

m = 5.1085 (ttT) . (3.4) 

Note that the supergravity analysis predicts not only the values of various correlation 
lengths, it also identifies which local operators in the dual field theory couple most directly 
to the lightest modes (in the strong coupling regime). In principle, it should be possible to 
confirm the predictions (3.2)~(3.4) via lattice simulations of the relevant correlators in AA = 4 
SYM — but this is not yet practically feasible. 



^''An open string does not source the two-form C at tree level, but fermion fluctuations on the string do 
couple to the RR field-strength. This implies that the e""*'^' tail in the C(x) correlator will have a coefiicient 
suppressed by 1/A relative to the short distance contributions. This shifts the center of the cross-over region 
by a tiny relative amount of order 0[(ln A)/A]. 
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4. Comparison to QCD 



Data for screening masses (i.e., inverse correlation lengths) in various symmetry channels, in 
SU{2) and SU{3) pure gauge theory, are available from 4d lattice simulations [28-31], as well 
as 3d simulations using the dimensionally reduced high temperature effective theory [32-36] . 
Results from both approaches are generally consistent with each other at the 10-15% level 
down to T = 2Tc. (For a good summary of the status of lattice calculations of thermal 
correlation lengths, see Ref. [37].) The dimensional reduction approach allows one to treat 
(with no additional computational complexity) QCD with any number of massless fermions. 
Using this approach, Ref. [32] reports results for many different symmetry channels and zero, 
two, three or four quark flavors. 

In Yang-Mills theory, the ratios of screening masses to the temperature arc very nearly 
temperature independent for l.bTc <T< AT, but (at least in the 0"^"'""'' and 0"* channels) 
decrease substantially as the temperature drops below 1.5Tc and approaches the confinement 
transition [28,29]. The degree of similarity between thermal QCD and A/" = 4 SYM is surely 
greatest within the window, l.STc < T < 4T, where screening masses scale linearly with 
temperature. Comparison of results for screening masses in SU{2) and S'f7(3) Yang-Mills 
theory reveal very little dependence on A'^c if A = g^Nc is held fixed. Consequently, the 
difference between iVg = 3 and A'^c = oo is expected to be quite small [32]. 

To make specific comparisons, we will use the results from Ref. [32] for T = 2Tc,^^ and 
we will also focus on the case of A/f = 2, where simulations were performed at a finer lattice 
spacing than for other (non-zero) values of A^f . QCD lattice simulations find that the smallest 
thermal screening mass is in the scalar channel, in agreement with the strong coupling 

AA = 4 SYM analysis. Table 2 shows results from Ref. [32] for the ratios of screening masses in 
different symmetry channels relative to the mass gap (the lightest screening mass) for Nf = 2 
QCD at T = 2Tc, together with the corresponding A = oo A/' = 4 SYM results from Table 1.^^ 
The same information is displayed graphically in Fig. 4. 

The lightest CT (or Rt) odd screening mass, which defines the Debye mass, is in the 
3%^* = Ol~ (axion) channel in both QCD and strongly coupled J\f = 4 SYM. As show in the 
above table and figure, there is essentially perfect agreement between the two theories for the 
value of the Debye mass to mass gap ratio, mD/wigap ~ 1.46. This remarkable agreement 
(to better than a percent) is surely somewhat fortuitous; in other channels the discrepancy 
between QCD and A/" = 4 SYM for these screening mass to mass gap ratios varies between a 
few percent and 30% (for the 2'*'+ channel). Overall, however, there is rather good agreement 
between the two theories. 

Actually, these results are for T = 2Ajjg. The ratio of Tc/Ajjg is close to one, but is not known with very 
high accuracy when A^f 7^ 0. 

^®In contrast to the (lighter) spin and spin 2 channels, Ref. [32] specified the time reflection symmetry but 
not the charge conjugation properties of the spin one operators they used. We have assumed that their spin 
one results should be compared with the lightest 3 = 1 channels with the given value of Rt which, for strongly 
coupled A/" = 4 SYM, are the spin one channels indicated in Table 2 and Fig. 4. 
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Table 2: Ratios of screening masses in indicated symmetry channels to the mass gap, for N{ = 2 
QCD at 2Tc, and A = oo, A/c = oo = 4 SYM. For both QCD and strongly coupled J\f=A SYM, the 
mass gap is the screening mass in the 0^'*' channel. 
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Figure 4: Ratios of screening masses in indicated symmetry channels to the mass gap (the 0++"'" 
screening mass), for N{ = 2 QCD at 2Tc (squares), and A = oo j\/' = 4 SYM (diamonds). The reported 
statistical errors on the QCD lattice results are smaller than size of the squares. 

Instead of focusing on ratios of screening masses, if one looks at their absolute size (in 
units of the temperature), then one finds that screening masses in QCD (at T ~ 2Tc) [32] are 
significantly smaller than in A = oo = 4 SYM: 



mgap _ ^ 


' 1.25(2), 
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= 2 QCD, T = 2T^; 




[2.34, 




= 4 SYM, A = oo, 


mo _ ^ 


'1.80(4), 


Ni 


= 2 QCD, T = 2T^; 




[3.40, 




= 4 SYM, A = oo. 



(4.1) 
(4.2) 

In A/" = 4 SYM, these masses (divided by T) are dimensionless functions depending only 
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on the coupling A. In the weak couphng regime, mga.p/T vanishes hnearly as A ^ 0, while 
itib/T scales as \/A. Presumably, both masses grow monotonically with increasing A and 
asymptote to the above A = oo, A/' = 4 SYM values. 

The smaller values for the Debye mass (and the mass gap) in QCD at 2Tc, as compared 
to A = cx) SYM values, naturally suggests that QCD plasma, in this temperature range, is 
most similar to AA = 4 SYM plasma at an intermediate value of A which is neither in the 
asymptotically weak, nor asymptotically strong coupling regimes. To test this hypothesis 
quantitatively, it will be necessary to compute subleading (in 1/A) corrections to A = oo 
screening masses in A/" = 4 SYM plasma. The Appendix discusses this in more detail, and 
reports some partial results on the first subleading correction to mo- 

Other possible extensions of this work include a comparison of correlation lengths as- 
sociated with flavor current correlators probing the dynamics of fundamental representation 
matter added to A/' = 4 SYM. In the gravitational dual, this corresponds to the addition of D7 
flavor branes as described in Ref . [38] . Lattice data is now available for a variety of mesonic 
correlation lengths in hot QCD [8,37,39,40]. We hope that this, and similar work, will 
clarify the degree to which A/" = 4 supersymmetric Yang- Mills plasma can mimic quantitative 
properties of the quark-gluon plasma produced in heavy ion collisions. 
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A. Prom weak to strong coupling 



The leading weak-coupling behavior of the Debye mass in A^ = 4 SYM is determined by 
the one-loop thermal correction to the gluon self-energy. One easily finds the lowest-order 
result [41,42] 

m^^^ = V2XT. (A.l) 

This is larger than the corresponding result for massless QCD (with A'c = 3) by a factor of 
V6 2.45 for Nf = 0, 3/a/2 ^ 2.12 for N{ = 2, or 2 for Nf = 3 (with Nf the number of quark 
flavors) . 

The next-to-leading order weak-coupling correction to mj) contains a logarithmic term 
whose coefficient may be easily extracted using perturbative effective field theory methods, 
together with a non-perturbative 0{XT) contribution which can only be computed via numer- 
ical lattice simulations. The analysis of Ref. [1] shows that changing the matter content of the 
theory from QCD to that of A/" = 4 SYM does not affect the next-to-leading order correction 
to the Debye mass (other than changing the value of m^^). The result is 

. X XT m^^'^ 

= + ^ 1^ !^ + ^ + 0{\'/^T), (A.2) 
47r Ai 

with the non-perturbative coefficient k = 0.64(2) for A^c = 3 [35].^^ The weak-coupling 
expansion of the Debye mass is much better behaved in A/" = 4 SYM than in QCD due to the 
larger value of . (The 0(AT) terms only become larger than m!^ when A exceeds 5.6.) 

The strong-coupling results for inverse correlation lengths derived in this paper have 
corrections which are suppressed by inverse powers of the 't Hooft coupling A, arising from 
a' corrections to type IIB supergravity. The first subleading corrections for the dilaton were 
investigated in Ref. [25] , with the result 

= 3.4041 [1 - 1.39 C(3) A'^/^ + . . . j (^3-) _ (^.3) 

To find the analogous subleading correction for the Debye mass, one cannot invoke the 
SL(2, R) symmetry of the tree-level IIB supergravity; an independent calculation of the axion 
mode is required. Corrections to the axion mass come from the original IIB action expanded 
around the a' corrected black hole background, as well as from new eight-derivative terms 
in the action affecting quadratic fluctuations around the original black hole background. 
Unfortunately these latter terms are not completely known at present. In principle these 
terms should follow by supersymmetry from the well known corrections, but the details 
are not yet fully understood. 

Possible corrections to the Debye mass arise from terms quadratic in the axion coupled 
to either the background curvature or the 5-form field strength. Ref. [43] has an explicit 

^''For SU{2), K = 0.63(2) [35]. Determinations of k from lattice simulations with Nc > 3 have, to our 
knowledge, not yet been performed, but the dependence of k on Nc is clearly very weak. Comparisons with 
four-dimcnsional simulations also show that, in pure Yang-Mills, 0{X^^^T) and higher contributions are not 
large compared to the 0(AT) terms when T > 2Tc [35]. 
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proposal for all the terms involving the axion and the curvature. The axion appears in the 
non-perturbative contributions to the modular functions /fc, but all those contributions are 
due to D-instantons and are exponentially suppressed by e~^/^^ ~ g-47riVc/A They are thus 
irrelevant in the large limit. The other two terms that are potentially important are 
12/0 R^DPDP and 6/2 R^{DPDP + DP DP), where P is proportional to the first derivative 
of the axion-dilaton pair. The axion derivative is the imaginary part of P. The factors /o and 
/2 are different functions of the axion-dilaton pair, but their leading weak coupling (i.e., large 
Nc) behavior is identical. In this limit, the axion contribution miraculously cancels between 
the above two terms from the proposal of Ref. [43] . Unfortunately the action proposed by 
Ref. [43] is not the unique SL{2, R) invariant extension of the known terms involving NSNS 
sector fields. In the absence of any rigorous arguments, an explicit calculation of the coefficient 
of the R^{d^a)^ term from a tree level 4-point function is needed. This computation has been 
performed recently in Ref. [44] with the result that the coefficient is non-zero, in contradiction 
with the proposal of Ref. [43]. However the precise normalization has not been determined, so 
a complete calculation of the sub-leading correction to the Debye mass is not possible at this 
time. In addition there will be non-vanishing terms involving the 5-form field strength, which 
is non-zero in the AdS black-hole background, and its derivatives. While some of these terms 
are known, see Ref. [45,46], there may be additional unknown terms involving the coupling 
of the 5-form field strength to the axion. 

We have, nevertheless, calculated the A"^/^ corrections that would follow from the simple 
conjecture of Ref. [43]. In that case the axion fluctuation satisfies the equation of motion 

5mx/=^5'"V'^5„a = 0, (A.4) 

where one has to use the a' corrected dilaton profile and metric in the lOd Einstein frame. 
Using the same shooting method as in Ref. [25], we find that the next-to-leading order (NLO) 
result for the Debye mass is 

mo = 3.4041 [1 - 0.52 C(3) A"^/^ + • • • ] (ttT) . (A.5) 

As with the dilaton, this correction leads to a smaller Debye mass as A decreases. Neglecting 
all further higher order corrections, this result becomes negative at A = 0.731. 

If one assumes, for the sake of discussion, that the additional A^'^/^ corrections due to 
eight-derivative terms in the corrected IIB action do not significantly change the result (A.5), 
then it is interesting to compare the strong and weak coupling expressions for the Debye 
mass. Figure 5 shows this comparison, along with one possible smooth interpolation between 
weak and strong coupling obtained by adding ncxt-to-next-to-lcading order (NNLO) terms 
to both Eq. (A. 2) and Eq. (A.5) and adjusting their coefficients to make the resulting curves 
tangent at a single point. 

1® Specifically, 0.6\^^^T was added to the weak coupling result (A. 2), and —5.0 A added inside the bracket 
of the strong coupling result (A.5). The resulting weak and strong coupling curves are tangent at A = 3.0. 
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Figure 5: Debye mass in A/" = 4 SYM as a function of A. The dashed (red) curve shows the NLO 
strong-coupling result (A. 5) based on the conjecture of Ref. [43] while the dotted (blue) curve shows 
the NLO weak-coupling result (A. 2). The solid (black) curve lying between these two NLO curves 
(for intermediate couplings) shows a smooth interpolation obtained by adding NNLO terms to both 
expansions, as explained in footnote 18, and adjusting their coefficients to produce results which are 
tangent at a single point. 

This choice of interpolation is certainly not unique, but it illustrates a reasonable choice 
given the available asymptotic information. This particular interpolation suggests that an 
AA = 4 SYM plasma will have a Debye mass which coincides with the T = 2Tc QCD result 
(4.2) when Asym ^ 2.5. The strong and weak coupling results would match better {i.e., 
smaller NNLO terms would be needed to make them interpolate smoothly) if the unknown 
additional contributions to the Debye mass from R^{d^a)^ terms in the effective action lead 
to a larger (more negative) correction at order A^'^^^. 
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